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Abstract 

We measured the low temperature subgap resistance of titanium 
nitride (superconductor, T c =4.6K)/highly doped silicon (degenerated 
semiconductor) SIN junctions, where I stands for the Schottky barrier. 
At low energies, the subgap conductance is enhanced due to coherent 
backscattering of the electrons towards the interface by disorder in the 
silicon (" reflectionless tunneling"). This Zero Bias Anomaly (ZBA) is 
destroyed by the temperature or the magnetic field above 250 mA 
or 0.04 T respectively. The overall differential resistance behavior (vs 
temperature and voltage) is compared to existing theories and values 
for the depairing rate and the barrier transmittance are extracted. 
Such an analysis leads us to introduce an effective temperature for the 
electrons and to discuss heat dissipation through the SIN interface. 

1 Introduction 

When a diffusive normal metal (or doped semi-conductor) (N) is electri- 
cally connected to a superconductor (S), electron-electron correlations (i.e. 
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a finite pair amplitude) are induced in the normal metal through Andreev 
reflections this is known as the proximity effect. The subgap conductance 
is a sensitive measurement of these correlations and a lot of recent works has 
been devoted to this subject, especially with the development of nanotech- 
nologies. In the original version of the proximity effect ||, the strength of the 
induced correlations depends on the barrier transmittance F at the interface 



and their extension is limited either by the thermal length Lt = <J 2 ^k B T or 



by Ly = y in a non-equilibrium situation (D is the diffusion constant 

of the normal metal); the phase breaking length L v or the geometry of the 
sample at small distances were not taken into account. 

Emergence of mesoscopic physics has enlightened the role of these effects, 
both when the contact between the superconductor and the metal is very 
good (high transmittance r«l) and in the opposite case (tunnel junctions). 
In high transmissive contacts, when Lt becomes larger than the typical sam- 
ple length L (supposed smaller than L v ), the re-entrance phenomenon takes 
place ||, which cancels out the standart proximity effect: instead of a steadily 
increase of the subgap conductance as T decreases, the conductance has a 
maximum at T such that Lt — L and decreases at lower temperature, 
theoretically recovering its normal value for T = OK (without interaction 
in the normal metal) [|]]. On the other hand, when the barrier is rather 
strong (low transmittance r< 1), the low temperature subgap conductance 
is small since the dominating Andreev reflection is a two-particle process 
which scales as T 2 . The subgap conductance for any value of < T < 1 
is described by the BTK theory ||, without taking into account coherence 
effects (ballistic normal metal). Coherent effects are also observed experi- 
mentally for r <C 1: when the Andreev reflected hole that traces back the 
incoming electron path is coherently backscattered towards the interface, the 
amplitudes for successive Andreev reflections add constructively. The sub- 
gap conductance is greatly enhanced as if the retroreflected hole did not feel 
the barrier ( "reflectionless tunneling" ) || . The conductance shows a peak at 
low energies, whose amplitude depends on the balance between the tunnel 
transparency and the rate of coherent backscattering. For a diffusive normal 
metal, the peak in the conductance increases when the coherent resistance 
increases and is maximum when the coherent normal resistance equals the 
barrier resistance. Eventually, if the coherent normal resistance increases 
above the barrier resistance, the situation evolves generically to the first case 
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(good NS contact) and the conductance shows a peak at finite energy 0. 

So far, reflectionless tunneling has essentially been observed in Super- 
conductor/Semiconductor junctions (S/Sm) (the only exception is the NS- 
SQUID devices by Pothier et al. [§]). In these systems, the low transparency 
of the interface is due to a Schottky barrier and to the mismatch of the 
Fermi velocities. The transparency of these junctions is intermediate (typi- 
cally Ts/sm — 10~ 3 ) between a very good interface (transmission coefficient 
T ~1) and an oxide barrier (r ~ 10~ 6 ). In semiconductors, annealing and 
surface cleaning processes or even the absorption of dopands by the electrode 
produce an increase of the sheet resistance of the semiconductor below the in- 
terface overlap. As reflectionless tunneling is a balance between the electron 
probabilities of crossing the barrier and backscattering to the interface, S/Sm 
systems are subject to this effect. But the superconducting material can also 
be weakened near the interface and consequently the BCS density of states of 
the superconductor can be smoothed by pair-breaking processes or creation 
of states below the gap. Zero-bias anomalies are always associated to such 
smooth conductance- volt age characteristics |7|, [|, [II], |IIf . Because this effect 
is very sensitive to the microscopic parameters near the interface as well as 
to the energy of carriers, one can take advantage of its observation both to 
feature the S/Sm contact and to investigate the thermalization of electrons 
in the normal part. The latter point has concentrated a lot of recent works in 
out-of-equilibrium mesoscopic normal conductors [0, [13|, [14], |15j . The prob- 
lem of carriers thermalization is even more crucial in S/N or S/N/S junctions 
because of the Andreev thermal resistance at the interface |L6], [17], [Tj|. It 
is also of practical importance for superconducting bolometers or Josephson 
field effect transistors. 

In this article, we report the observation of zero-bias anomalies in tita- 
nium nitride/heavily doped silicon junctions at very low temperature (down 
to 30mK). In a first part, the samples are characterized with various measure- 
ments (contact resistance, sheet resistance, weak localization). In a second 
part, the differential resistance of the junctions is measured as a function of 
temperature, voltage and magnetic field. The third part is devoted to the 
quantitative comparison between the observed zero bias anomaly and the 
theory for a planar SIN junction. A good agreement is obtained for tem- 
perature behavior, allowing us to extract the depairing rate and the barrier 
transmittance of the junctions. However, discussion on the voltage response 
of the junctions leads us to consider the effective temperature of the carri- 
ers in the silicon, which is well above the phonon bath temperature. This 
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overheating effect is discussed in the fourth part. 



2 Samples and materials characterisation 

The samples are fabricated from a TiN(1000A)/Si n ++ (P doped on 
dsi:p=0.6fim) bilayer deposited on a Si 8' wafer. First, the wafer is oxi- 
dized on 13nm. Phosphorus is then implanted at 15 keV, 2.10 15 cm -2 , fol- 
lowed by a recristallisation heat treatment at 650° C and a 30 minutes activa- 
tion/diffusion treatment at 1050°C with oxygen. The substrate is deoxidized 
and a lOnm Ti and lOOnm TiN bilayer is deposited. Optical lithography 
defines a Tranverse Length Method (TLM) pattern with various distances 
L — 1,2, 5, 10, 20, 50, 100, 200 and 500/zm between large TiN pads (typically 
lOOOx 1000/zm 2 ) (see bottom insert of figure [IJ). The TiN/Ti is etched and a 
final heat treatment at 720°C in N2 atmosphere during 20 seconds provides 
TiN densification and forms a 40nm thick TiSi2 layer. Nevertheless, the 
zero bias anomaly reported in this work appears rather insensitive to this 
heat treatment. On some samples, another etching is performed to define 
two TiN pads fingers with lateral dimensions w=10/xm facing each other and 
connected to larger TiN reservoirs. 

Although titanium nitride has been used for many years in microelec- 
tronics as a diffusion barrier, ohmic contact and gate electrode in field effect 
transistors J0|, its superconducting properties have been thoroughly studied 
only recently f20| . Its transition temperature is T c =4.6 K for our samples. 



STM measurements give a gap of 250/iV at T=1.4K [21]. The relation be- 



tween the superconducting gap and the transition temperature departs from 
the BCS theory probably because titanium nitride is a granular supercon- 
ductor. Its room temperature resistivity is 85/zf2.cm. 

Doped Si:P has been studied for a long time because of its great impor- 



tance in microelectronics. Alexander et al. |22| evaluated the Mott-transition 
donor concentration n c = 3.10 18 cm~ 3 and the concentration at which the 
Fermi level of the electron system passes into the conduction band of the 
host crystal = 2.10 19 cm -3 . Heslinga et al. [^3] studied the inelastic 
lifetime in heavily doped Si:P. They found l/r in (s) = 1.1 10 9 T(K) 2 - 2 at 
n = 2.10 19 cm -3 from T=1.2 to 4K. In our samples, doped silicon Si:P forms 
the normal part of the junction with a donor concentration n e =2.10 19 cm -3 
over a depth of P implantation dsr.p=0-Q^m. From the TLM geometry, we 
estimate the sheet resistance of the doped silicon and the resistance Rnn of 



4 



the interfaces at T=4K: 

R = ^ + 2R NN (1) 

where i? n is the sheet resistance of the doped silicon, N n = w/L the number 
of squares of silicon in parallel between the two SIN contacts, L' and w 
the length and width of the TiN overlap and L the distance between the 
two TiN electrodes (see inset figure |l|). Following Giaever's calculation f24| , 
Rnn = \/RbRh/ w is the normal state resistance per contact (determined 
at temperatures just below T c , to eliminate the resistance of the TiN pads), 
with Rb the barrier resistance expressed in Q./im 2 and Rq the sheet normal 
resistance below the overlap (eventually R^ > Rn for the bare film) . From 
equation [1], we see that the resistance is linear with the distance L. This is 
what we observed at small L on figure |T[ For L > 50/im, dispersion of the 
current lines on the sides of the normal part of the S/N/S system becomes 
noticeable (the TiN pads are on top of an infinite silicon layer). Therefore, 
we can recover the experimental curve by adding in this case two squares of 
silicon to N n . 

The linear fit at small lengths (using w = 1000/zm) gives at 4K (see inset 
figure |l|) : 

R(tt) = 0.024 L(fim) + 0.156 (2) 

We deduced a sheet resistance of the bulk silicon (between the interfaces): 
R n = 24f2. The resistivity at 4K is then p=14.4 /iQ.m in very good agree- 
ment with reference ||22|| . The normal resistance of each interface is Rnn 
7.x 10- 2 Vl for w=1000/im. 

Given the effective masses of silicon m* = 0.321m e (with m* t = 0.19m e 
and m\ = 0.916m e ) and the valley degeneracy N=6, we used a free electrons 
model to calculate the parameters of the doped bare silicon: 

k F = [-Y^J = 4 - 62 iO 8 ™ -1 and X f = l3.Qnm (3) 
hkp 

C e = — - — = Q.Qnm and kp£ e = 3 (4) 
pe 2 n e 

D = -v F £ e = 3.67 10- 4 m 2 .s _1 (5) 
3 

Those values ensure that the doped silicon is in the metallic regime, since 
kp£ e > 1 and n e > n c . 

Finally, we measured the magnetoresistance of the bulk silicon using a 
long and wide bar (20mm x 1mm) of Si n ++ , at various temperatures (figure 
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Figure 1: Tranverse Length Method (TLM) measurements at 4K for different 
lengths L. In upper inset, the same curve for small L. We deduce the sheet 
resistance of the bulk silicon R D = 24fi and the normal state resistance of 
an interface Rnn = 7.8 lCT 2 f2 for a w=1000/im width junction. The dashed 
curve is obtained by adding two squares of silicon in parallel to take into 
account the conduction through the silicon layer outside the gap between 
the TiN electrodes. In bottom inset, schematic view of the samples. 
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0). We fitted the experimental curves with the 2D and 3D theories of weak 
localization 112511 : 



O — ^Boltzmann + A<7 (6) 

Aa = ^/ 2 (2^f ) at 2D (L v > d Si:P ) (7) 
-/ 3 (2^) at 3D (L v < d 5i:P ) (8) 



with 



/ 2 (a;) = ln(x) + 



n=0 



Lh = yfrj (2eH) and are the magnetic and the phase-breaking lengths 
and dsi-.p the depth of doped silicon. ip(x) is the digamma function. We 
deduced t v = L^/D for various temperatures (see table |l|). As we were 
not able to fit the experimental curve at T=700mK, we concluded that the 
cross-over between the 2D and 3D regimes lies in this temperature range 
(corresponding to L v ~ dsi-p). For lower temperatures, we used the 2D 
theory and for T=1.4K, the 3D theory. At T=1.4K, we measured t v = 
0.48ns, in good agreement with Heslinga et al [23] who obtained t v = 0.58ns 



at T=1.2K. At lower temperature, these results can be compared to the 



expression given by Altshuler et al. [ 26| , where only R n enters as a parameter: 



which gives t v = 1.3ns at IK (see table |l|), in relative good agreement with 
our experimental results. 

As distances between the TiN electrodes are larger than or of the order of 
the phase-breaking length for all studied samples, as far as the coherence is 
concerned, one can divide the samples into three parts: the silicon resistance 
and the two interfaces. No coherence effect should link both interfaces and 
they could be studied as two N/S systems in series. This is confirmed in the 
next section. 
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Figure 2: Magnetoresistance of a wide wire of Si n ++ measured for various 
temperatures. We deduced 1^=420, 620 and lOOOnm respectively at 1400, 
360 and 240mK. In inset, the magnetoresistance at 700mK shows the cross- 
over between 2D and 3D regimes (curves at L v ,=<isj : p=600nm depth of doped 
silicon). 
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Temperature (mK) 


r (ns) (exp.) 


(ns) (calculated eq. ||) 


1400 


0.48 (3D) 


0.9 


700 


cross-over 2D-3D 


1.9 


360 


1 (2D) 


3.6 


240 


2.7 (2D) 


5.4 



Table 1: Coherence time for various temperatures deduced from figure |2| and 
reference [26]. 



3 Zero bias anomaly in TiN/Si planar junc- 
tions 

We measured the sample resistance versus temperature (down to 30mK) 
by a standard four probes lock-in technique (I ac = 10 nA at 180Hz). We 
studied samples of different lengths (from L=lyum to 500/iin) and they all 
show the same behaviour at low temperature (see inset figure |3|). Just above 
4K, we observe a step in the resistance due to the superconducting transition 
of the TiN electrodes. Independently, by measurements on a long and wide 
bar of doped silicon, we have checked the silicon resistance does not depend 
(or only slightly) on temperature and voltage (see figure |2] and part ||) . Since 
we know the normal resistance from TLM measurements, we can plot the 
resistance per contact: R C (T,V) = ^(Rtotai(T,V) — Rsi) where Rsi is the 
resistance of the silicon between the two contacts. So we are able to plot the 
resistance per contact R c as a function of the voltage drop at the interface 
Vi. The voltage drop at each interface is Vi — \{Vtotai — RsJ) where / is the 
applied DC current. 

In the rest of the paper and in figure we plot R c . Below T=4.2K, 
the resistance increases as expected for an SIN junction. Anticipating the 
discussion in the section 4, its temperature dependence is well fitted within 
the BTK model || above T ~ 400mi^. The resistance increases only by a 
factor of 6 between 4K and 400mK, indicating a strong departure from the 
sharp BCS density of states. This is taken into account via a high damping 
factor Tg in the superconductor (Ts/A = 0.11 in the BTK adjustment) 
while the barrier is rather opaque (transparency 3.4 10 -2 ). At 250 mK, 
it shows a maximum and then decreases. This can not be explained as a 
precursor of a Josephson current, since this effect is observed for distances 
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Figure 3: Resistance per contact versus temperature at zero bias for a 20 x 10 
fim 2 sample. The fitted curve used the following parameters. BTK (dashed 
curve): A = 0.21meV, T s /A = 0.11, Z=5.3 and R NN = 7.45ft. Volkov 
(solid curve): A = 0.22meV, T s /A = 0.15, 7/A = 0.27, e b /A = 0.009 
and Rnn = 7.45ft. In inset, total resistance-temperature characteristics for 
samples of various lengths L and w = 1000/im: from top to bottom, L=l, 
200 and 500/im. 



between superconducting electrodes much larger than the coherence length 
L T = v /hD/(2nk B T) = 120nm at 30mK (see figure D). We interpret this 
decrease of the resistance at low temperature as due to reflectionless tunneling 
1,1- 

We also measured the differential resistance versus DC voltage at different 
temperatures (see figure ^) by applying a DC and ac current and by recording 
both the DC and ac voltages. 

Except at low voltages and low temperatures, the R C (V) characteristics 
show the same smooth behavior for all samples: as the voltage is increased, 
the resistance steadily decreases and reaches its normal value above the gap 
voltage: there is no shoulder around eV = A. No shoulder is observed either 
at the gap voltage for previously studied Sm/S contacts exhibiting reflec- 
tionless tunneling Q. At low temperature, the resistance shows a dip (ZBA) 
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consistent with the reflectionless tunneling regime. When the temperature 
is increased (see figure (D, the ZBA amplitude vanishes and disappears at 
250mK. Moreover, as a function of the voltage, the ZBA disappears on the 
same energy scale than as a function of temperature, i.e. roughly 20/iV. 

Finally, we measured the differential resistance versus the DC voltage for 
various magnetic fields (see figure |5|). Two behaviors can be distinguished. 
When the magnetic field is perpendicular to the S/N interface (along y, 
see figure [3]), the resistance decreases at very small field and the ZBA is 
destroyed at 30G. On the contrary, for weak applied magnetic field parallel 
to the interface (along z), the resistance background is unchanged at small 
field and only the ZBA diminished. Then, above 200G, the overall resistance 
decreases. The ZBA is divided by a factor two for H* ~ 200G and completely 
disappears at 400G. 

4 Proximity effect in planar SIN junctions 

Many theoretical works have been completed since Kastalsky experiments 
p7| , p8| , p9| . Volkov |30] considered planar SIN junctions and, using a the- 
ory based on the Usadel equation, he calculated the conductance-voltage 
characteristics of this system at any voltages and temperatures. Assuming 
e b << 7 in , with e h the barrier energy (see below) and 7j n the depairing rate, 
he obtains in the case of long junctions (L » L b ): 




J de\A{e)f N , z {e,y) 



(10) 




with 



\A{e) = Re- 



e b %n A 2 



a*(A-|e|) 



2 + 7^A 2 -( e + T\) 




(11) 



h 



(12) 



Tin 



in 




(13) 
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Figure 4: TOP: Resistance per contact versus voltage for different tempera- 
tures for a 20x10 fim 2 sample and the theoretical curve at 30mK after Volkov 
model (parameters extracted from the temperature dependence of R(V=0)). 
BOTTOM: Resistance per contact versus voltage at large voltages, BTK fit 
with T=320mK (dotted line) and BTK fit with a Wiedemann-Franz law (full 
line) for a 1x10 fim 2 sample. The BTK parameters are also extracted from 
the temperature dependence of R(V=0): A = 0.22meV, r s /A=0.12, Z=5.3, 
R NN = 7.45fi and T = 320mK. 
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Figure 5: Resistance per contact versus voltage for various magnetic fields 
applied parallel (top figure) and perpendicularly (bottom figure) to the in- 
terface for a 1x10 fim 2 sample. In inset of the top figure, the depen- 
dence of the pair-breaking energy versus the applied magnetic field is 
obtained from the evolution of R(V=0) versus H. The fit gives the pair 
breaking parameter 7 deduced from the evolution of R(V=0) versus H: 
7i„(#)/A = 0.405 + 4.5 10~ 7 H(G) 2 §(J. 
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with A the superconducting g small damping in the superconduc- 

tor, Ti n the inelastic time and the barrier length, Rb the tunnel bar- 
rier resistance and Rn the sheet resistance of the normal part. f^, z = 
i[tanh(|^0 — tanh(|^0] is the equilibrium distribution function. 

The first term in \A(e) (equation |TT]) describes the pair amplitude in the 
normal metal. It is non-zero (although small) in spite of the null value of the 
electron-electron interaction potential in this region: this is the proximity 
effect. As e& — > (the barrier is more and more opaque), the proximity effect 
reduces and the pair amplitude vanishes. In the same way, this pair amplitude 
is zero when there is no coherence effect (7^ — > +00). The second term is the 
BCS normalized superconducting density of states, which leads to the usual 
conductance- voltage characteristics (tunnel Hamiltonian approach). At zero 
temperature and voltage, Gns = Gnn\J '^b/lin as long as e^, 7 in << A. The 
width of the peak is roughly proportional to j in . 

As we said in the introduction, S/Sm junctions are intermediate between 
good metallic junctions and oxyde tunnel barriers. We choose to treat our 
junctions within the SIN Volkov's model, taking into account that the abso- 
lute contact resistances are high and that estimation for the Schottky barrier 
gives poor interface transparencies. We will show that applying Volkov's 
model forces us to attribute a much larger sheet resistance for the silicon 
layer below the overlap (R^=315 Q) than for the bare silicon (i? n =24 Q). 
The depairing rate below the interface is also found to be much lower than 
the phase breaking rate for the bare film. Within this theory, we stress that it 
is impossible to consistently describe the ZBA taking the silicon parameters 
in the bare film unchanged below the overlap. This illustrates how the full 
quantitative analysis of the ZBA provides distinctive informations about the 
microscopic parameters of the S/Sm contacts, otherwise unaccessible. 

Figure |^ shows the experimental data for the zero bias resistance versus 
temperature and the best Volkov's model adjustment. The parameter are 
the following: titanium nitride superconducting gap A=0.21meV, damping 
factor rs/A=0.12, normal resistance Rnn=7-4:5Q, barrier transparency per 
channel T = 3.4 10~ 2 , and finally the pair breaking energy 7j n = 60 — 85fieV. 

There are two separate sets of parameters in this list: A, T s , Rnn are 
determined in the high temperature range (and also independently by the 
large voltage range in the differential resistance curves, see later on); -y in and 
T are fixed by the low temperature range of the curve. In principle, V enters 
also as a parameter in the high temperature range, but due to the large value 
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of r§, any relative small value of T satisfies the adjustement. 

In the high temperature range (above 500mK), we also fit the resistance 
versus temperature curves with the BTK model || , using the same values for 
the parameters than for the Volkov's adjustment (see figure |3|). Because in 
this range of energy the phase coherence is likely negligible, we find that both 
formalisms describe equally the experimental data. Below 500mK, coherence 
effects appear and the BTK model is no longer valid. 

We now comment on the various parameters: the damping factor Ts is 
large probably because the disorder weakens the superconductivity at the 
interface. In particular, TEM image of the contacts indicates that a TiSi2 



40nm thick layer is produced by the thermal treatment at the interface |3T 
From the transparency, we estimate the barrier resistance Rj, = A ^ Xp ^ ~ 
18Q. fim 2 . From the absolute value of Rnn, we can deduce the sheet resis- 
tance of the silicon underneath the interface: R% = Rn £ w =3150. This value 
is 15 times larger than the bulk value obtained with the TLM measurements 
(i? n =24 Q). Note that without the determination of T by the Volkov's fit of 
the ZBA, another set of T — R^ parameters could explain the absolute value 
of R NN , namely for instance R* n = R D = 2AQ and T = £ { -^^- = 2.35 1(T 3 
(with Xp = 13.6nm). This set will give a good BTK fit above T = 500m.fr, 
but the T is too small to give enough Andreev reflection to provide the ZBA. 
The ZBA is explained only if both the normal sheet resistance of the silicon 
layer under the overlap and the transparency of the Schottky barrier are large 
enough. As said previously and also in reference 0, the semiconductor layer 
under the overlap could have a much larger sheet resistance than the native 
layer and is close to the metal/insulator transition (kp£ e — 1)- In such a 
case, the phase breaking length should also strongly decrease in the vicinity 
of the overlap as compared to the bare film. 

Moreover, we obtain = 60 — 85fieV, which leads to r in = h/^ in =9>- 
llps. This value, which is an upper limit for t v , is smaller by two orders of 
magnitude as compared to t v ~2. 7ns (L^=l/im at 240mK) estimated from 
magnetoresistance measurements in the bulk silicon layer. Interestingly, such 
small Ti n has been reported in tunnel experiments in copper wires || or in 



gold wires from analysing the Josephson effect in SNS junctions ||13|| , and 
the discrepancy with phase breaking times measured by weak localization 
measurements also quoted. We do not know if it is coincidental, or if local 
measurements of r, n near an interface generally lead to such underestimations 
as compared to bulk measurements. 
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In summary, the precise analysis of the reflectionless tunneling effect leads 
us to deduce that the ZBA is not due to coherent backscattering over long 
distances (about L^—lyiva) into the bulk silicon, but to coherent backscat- 
tering in a short disordered layer. 

Finally, we discuss the effect of an applied magnetic field parallel or per- 
pendicular to the junction. First, we measured the resistance-voltage charac- 
teristics at 30mK for various magnetic fields parallel to the interface (along 
z, see figure [I]). We note that the zero bias anomaly is divided by a factor 
two for an applied fields H* ~ 200G and completely disappears at 400G. 
We can estimate the field necessary to put one flux quantum $ = h/e in a 
square of size L\ n . Following Marmorkos et al p2| , it gives the critical field 
to destroy the zero bias anomaly. Under the interface, we deduce from the 
fit of our experimental curves, r in =8-llps, so L in = y/ Z}r irt =54-64nni with 
D = 3.67 10~ 4 m 2 .s _1 . Consequently, H c = <& /L 2 n ~1T, which is much 
higher than the observed values of the field. Volkov et al. |33j proposed 
another mechanism to destroy the coherence of the Andreev pairs. The mag- 
netic field leads to a screening current, which gives a dependence of the phase 
of the order parameter with the coordinate z. This situation is comparable 
to Andreev interferometers ||: if the two electrons encounter a supercon- 
ducting phase varying over ir, destructive interferences occur and the ZBA 
disappears. Then, the depairing depends on the magnetic field through: 

j in (H)/A = 7i n (0)/A + HD(nH\/$ ) 2 (14) 

with A the London penetration depth of the superconductor. Using Volkov 



theory (see equation [10]), we calculate the value of the depairing rate at 
zero voltage for various fields and obtain the fitting curve: ji n (H)/A = 
0.405 + 4.5 10~ 7 H(G) 2 (see inset figure §), which leads to A = 790nm. This 
value seems rather high compared to A ~ 200nm for NbN. To calculate this 
penetration depth, we take the diffusion coefficient in the bulk silicon (see 
equation [5]) which maybe lower under the interface because of the disorder in- 
duced during annealing. Then, the penetration depth maybe overestimated. 
Nevertheless, this mechanism may cause the destruction of the interferences 
in our sample, since they are not sensitive to the decoherence induced by flux 
in electron-hole trajectories in the normal part (critical field too high). 

But great care should be taken with the orientation of the magnetic field. 
We also measured the resistance-voltage characteristics at 30mK for various 
magnetic fields perpendicular to the interface (along y). The overall subgap 
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resistance decreases for very small fields. Since titanium nitride is a type 
II superconductor, we attributed this effect to the appearance of vortices: 
the total resistance decreases because the junction normal resistance due to 
the vortices is less than the superconducting junction resistance. Secondly, 
the zero bias anomaly is divided by a factor two for a weak applied field 
H% ~ 30G (see figure ||). We can estimate the demagnetization factor under 



the interface. According to Zeldov et al. |[34|| , the field under the interface 



is given by H interface ~ y 'w / ' d TiN H appUed ~ 10H appUed , with d TiN = lOOnm 
thickness of the TiN film and w = lfim. This value agrees with the ratio of 
characteristic field: H^/H* ~ 7. Consequently, the ZBA for magnetic field 
applied perpendicularly to the interface disappears for weak field because of 
the demagnetization factor of the superconducting film. 

The vortices may also explain the decreasing around 200G of the resis- 
tance background when the magnetic field is parallel to the interface: if the 
field is not strictly parallel to the interface, some vortices may appear with 
a small perpendicular component of the field. 



5 Electron heating and effective temperature 

The Volkov's theory or the BTK theory are not able to fit our experimen- 
tal resistance- voltage curves if we suppose an equilibrium Fermi distribution 
with a base temperature T . Paradoxically, we note that the maximum of 
differential resistance happens precisely at a voltage Vo — 20/iV such that 
eVo — ksTo, where To ~ 250mK is the temperature at which the zero bias 
resistance is maximum. Such coincidence happens also in the context of the 
re-entrance effect |3], |35|] and has been observed in most of the reflect ionless 
tunneling experiments [|7|, P| fTI]. But this is not predicted in the theoretical 
models. From the model, we expect that the voltage brings to a much smaller 
suppression of the ZBA as compared to experimental data (see figure |4j). A 
way to restore a good accordance between the model and the experimental 
voltage characteristics is to choose the effective temperature of carriers as 
an adjustable parameter [3^J. At a fixed voltage, we find the elevated tem- 



perature such that the model gives precisely the measured resistance value. 
Therefore, we construct the variation of the effective temperature T e ff as a 
function of the applied voltage (figure |]). We attribute this temperature to 
the carriers in the doped silicon film in between the two TiN electrodes. 
What is happening first at very low temperature when a finite voltage 
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is applied to the TiN/Si n ++ /TiN sample? Some Joule power is dissipated 
in the silicon part which is very hardly evacuated either in the phonon bath 
or in the contacts because the electron-phonon coupling is very small at low 
temperature and the Andreev thermal resistance at each N-S interface is 
very large Q. Consequently the electrons inside the silicon are overheated. 
Overheating has two causes. First, the Andreev thermal resistance, which 
depends exponentially on the effective temperature, is mainly responsible for 
the elevated temperature at low voltage, and one can neglect the gradient of 
temperature in the silicon between the superconducting contacts. Secondly, 
at high voltage, this gradient is non negligible and given by the Wiedemann- 
Franz law. We suppose that the inelastic electron-electron scattering time 
is short as compared to the length of the sample (1 — 20/im), to have an 
quasi-equilibrium Fermi distribution with an effective temperature. 

At low voltage, below 20 /xV, the electronic temperature increases very 



rapidly up to 320mK. Hoss et al. [IS] observed a very similar behavior in 



lfMm long Nb/Au/Nb samples at low currents. To understand this results, 



we used a BTK-based model of dissipated power fL8| . The power dissipated 
through a NS interface is given by: 

P(T e , 7) = — / de e [/(-r^) - /(tt^)][1 - A(e) - B{e)\ (15) 

e J-oo KBJ-e 

with T e and T electronic and phonon temperatures, Vi voltage drop at the 
interface, A and B the Andreev and normal reflexion probabilities [0, Gnn — 
1/Rnn the normal conductance. We fit the curve with the parameters used 
earlier in the BTK fit: A = 0.21meV, T s /A = 0.14, Z = 5.3 (which gives a 
barrier transparency of 3.4 10~ 2 ), Rnn = 7.45f2 and T = 30m.fr. We assume 
that all the electric power is dissipated through the two NS interfaces, i.e. 
that the electron-phonon length at this temperature exceeds the distance 
between the two superconducting interfaces. Then, by equating P(T e ,Vi) = 
Vtotail, with Vtotai the total voltage and / the current across the sample, we 
obtain T d (Vi). 

These parameters describe the experimental points of figure |6| at low 
voltages. Since the thermal resistance of the N/S interface decreases ex- 
ponentially with temperature, for larger effective temperatures the heat is 
rapidly evacuated in the superconducting electrodes and a thermal equilib- 
rium is reached, giving the saturation of the effective electronic temperature 
supposed constant in the silicon to around 300mK (solid symbols in figure 
©• 
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Figure 6: Electronic temperature (□ points) deduced from comparison be- 
tween theoretical (Volkov's model with previous parameters for R(V) curves; 
see text) and experimental R(V) curves of the 1x10 /im 2 (top) and 20x10 
jum 2 (bottom) samples. The solid curve is the Wiedemann-Franz law 

T e (V) = y/Tj. + |^(g) 2 with T Si = 320mK. The diamond-curve is de- 
duced from the dissipated power model, with A = 0.2lmeV, T s / A = 0.14, 
Z = 5.3, R NN = 7A5tt and T = 30mK. 
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At higher voltages, the effective temperature continues to increase. Our 
hypothesis of a constant effective temperature in the silicon should fail. 
The simplest analysis is to consider the Wiedemann-Franz law, globally 
for the whole TiN/Si n ++ /TiN sample. In the 1x10 fim 2 sample, it fol- 
lows a Wiedemann-Franz law (interacting hot electron regime) T e (V) = 



T"si + i^(f^) 2 ) w ^h Tgi = 320mK temperature of the electron at the 
interface and V the total voltage applied to the SNS system (see figure H). 
In the 20x10 /j,m 2 , the effective temperature is less than predicted by the 
Wiedemann-Franz temperature. Heslinga et al [23| give an approximation 
of L e _ p h ~ 2fim at 320mK. This value is intermediate between the length of 
our two samples and explains their different behaviors. In the 20x10 fim 2 
sample, as the length of the sample matches the electron-phonon length, 
electrons can interact with phonons and be cooled under the Wiedemann- 
Franz temperature, whereas in the 1x10 /im 2 , electrons are only heated by 
electron-electron interaction and the Wiedemann-Franz law is valid. 

Finally, we test our estimation of the effective temperature using the BTK 
model. On figure |4] (bottom), we observe a much better accordance with the 
data if we introduce the effective temperature given by the Wiedemann- 
Franz law than if we use the base phonon temperature. For consistency, all 
the parameters except the temperature are obtained from the BTK fit of 
R(V = 0) versus temperature. 



6 Conclusion 

In conclusion, we measured the differential resistance versus tempera- 
ture, applied voltage and magnetic field of semiconductor/superconductor 
junctions. TiN/Si n ++ heterostructures proved to be a new and interest- 
ing tool for the study of proximity effect. We observed a zero-bias anomaly 
due to reflectionless tunneling. By comparing our results to a proximity ef- 
fect theory |30| , we find that the sheet resistance of the silicon underneath 



the interface is much larger than the bulk value and that the quasiparticule 
lifetime is much shorter than in bulk silicon. We explain these discrepan- 
cies by disorder induced by annealing during the process. Interestingly, the 
quasiparticule lifetime is comparable to the electron-electron interaction time 
deduced from tunnel experiments in copper wires |§ or in gold wires from 



analysing the Josephson effect in SNS junctions At very low tempera 
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ture and finite voltage, the effective electronic temperature in the silicon is 
much higher than the phonon temperature. This is valid for separation be- 
tween superconducting contacts as large as 20/xm. The effective temperature 
increases very rapidly at low voltage and then follows a Wiedemann-Franz 
law. This behavior is well explained by a model of heat dissipation through 
a N/S interface. 

This rapid raise of the temperature of the carriers with the injected power 
is one of the interesting properties of this system. It could be used to make 
a bolometer, by measuring directly the zero bias conductance as function of 
the absorbed power in the silicon J?6| . 
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